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I. INTRODUCTION 

The gluon self-coupling in QCD opens the possibility of 
existing bound states of pure gauge fields known as glue- 
balls. Glueballs (G) are predicted by several theoretical 
formalisms and by lattice calculations. For a comprehen- 
sive review on the current status of theoretical and ex- 
perimental aspects of glueball studies we quote Ref. [l| 
and Ref. [2j , respectively. Many mesons have stood up as 
good candidates for the lightest glueball in the spectrum 
and in particular the scalar sector (J PC = ++ ) seems 
promising. The mesons /o(1500) and the /o(1710) have 
been considered the principal candidates for the scalar 
glueball [3, Q. However, in this mass region the glue- 
ball state will mix strongly with nearby qq states [J, [5J . 
More recently, the BES collaboration observed a new res- 
onance called X(1835) |6j. It is an important candidate 
for glueball and the nature of meson X(1835) has several 
interpretations. One of them consider it a pseudo-scalar 
glueball (J PC — _+ ) as first suggested in Ref. [7| and 
afterwards in 8] . 

Recently, the clean topologies of exclusive particle pro- 
duction in electromagnetic interactions hadron-hadron 
and nucleus-nucleus collisions mediated by colorless ex- 
changes such the QCD Pomeron or two photons have 
attracted an increasing interest Q. The cross sections 
for these processes are smaller than the correspondent 
inclusive production channels, which it is compensated 
by a more favorable signal/background relation. Experi- 
mentally, exclusive events are identified by large rapidity 
gaps on both sides of the produced central system and 
the survival of both initial state particles scattered at 
very forward angles with respect to the beam. 

Here, we will focus on exclusive glueball production 
in two-photon and Pomeron-Pomeron interactions in co- 
herent nucleus-nucleus collisions at high energy colliders 
(RHIC and LHC). In these cases, the photon flux scales 
as the square charge of the beam, Z 2 , and then the cor- 
responding cross section is highly enhanced by a factor 
a Z 4 k 10 7 for gold or lead nuclei. A competing chan- 
nel, which produces similar final state configuration, is 
the central diffraction (CD) process. Such a reaction is 
modeled in general by two-Pomeron interaction. Exper- 



imentally, the separation of these channels is somewhat 
difficult and from theoretical point of view the Pomeron- 
Pomeron are subject to large uncertainties at collider en- 
ergies. One goal of present work is to compare the cross 
sections for these two channels in the production of glue- 
ball candidates. This paper is organized as follows: in 
next section we present the main expressions for cross 
section calculation of two-photon and Pomeron-Pomeron 
processes and in last section we shown the numerical re- 
sults and discussions. 



II. CROSS SECTION CALCULATION 

Let us start with the glueball production in photon- 
photon scattering at coherent heavy ion collisions using 
the Weizsacker - Williams approximation (EPA approxi- 
mation). In such an approach, the cross section for a two 
quasi-real photon process to produce a glueball state, G, 
at center-of-mass energy W 77 factorises into the product 
of the elementary cross section for 77 — > G convoluted 
with the equivalent photon spectra from the colliding ions 
91: 



ct 77 (>L4->- A + G + A) 



dki dki dn-y dn^ 
k\ ki dk\ dki 



<r( 77 ^G),(l) 



where k\^ are the photon energies and dn/dk is the 
photon flux at the energy k emmited by the hadron A. 
The photon energies determine the center-of-mass energy 
W 77 = \J\k\ki and the rapidity Y of the produced sys- 
tem. Namely, one has fci,2 = (W77/2) exp(±F) and Y — 
(1/2) In (fei/fe). In addition, ^snn is the center-of-mass 
energy of the ion-ion system and the Lorentz relativistic 
factor is given by 7l = • v /sj\r./v/(2m./v)- In particular, in 
the numerical calculations we use ^snn = 0.2 (5.5) TeV 
and 7l = 109 (2930) for RHIC (LHC). 

In the EPA approximation, the flux of equivalent pho- 
tons from a relativistic particle of charge Z is determined 
from the Fourier transform of its electromagnetic field. 
For an extended charge with electromagnetic form fac- 
tor, Fa(Q 2 ), the energy spectrum can be computed as, 



dn y/A (x) _ aZ 2 A(x) f Q 2 -Q 



dk 
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where x — k/E is the fraction of the beam energy carried 
by the photon and A(x) = 1 — x+(l/2x 2 ). Moreover, a — 
1/137 and Q 2 is the four- momentum transfer squared 
from the charge, with <3 2 nin ~ (a:mjv) 2 /(l — x). 

The glueball production in two-photon fusion can be 
calculated using the narrow resonance approximation 



10]: 
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where T(G — > 77) is the partial two-photon decay width 
of G, Mq is the glueball mass and J is the spin of the 
state G. Here, we compute the production rates for 
the mesons / (1500), /o(1710) and X(1835) [ll|, respec- 
tively. The reason is due to they have been mentioned as 
possible glueball candidates by phenomenologists [l|, [2| ■ 

Some important comments are in order. The predic- 
tions for the two-photon component are in practice some- 
what difficult as the branching ratios have not been mea- 
sured. To compute numerical values for the meson (glue- 
ball) cross section in two-photon reactions estimates for 
the two-photon decay widths are needed. The determi- 
nation of them depend upon whether the meson state is a 
pure quarkonium, pure gluonic or a mixed hybrid state. 
For a pure quarkonium state the width can be related 
(at leading order) to the two-gluon width, T(qq — > gg). 
Namely, T(qq ->• 77) ~ D c e\ (a/a s ) 2 T(qq -t gg), where 
D c = 9/2 is the colour factor and e q is the relevant quark 
charge. One can estimate the two-gluon width from the 
total width for the meson state and the theoretical expec- 
tation that the qq — > gg branching ratio 1 is of order a 2 
[l2| ■ In case of a pure gluonic state, the two-photon width 
can be computed using a nonrelativistic gluon bound- 
state model as performed for instance in Ref. [13| . There, 
the unknown parameters as the digluon wavefunction, or 
its first/second derivative at the origin, are determined 
by using measured values of T(J/tp — > G7). 

Now, we compute estimates for the two-photon widths 
assuming pure qq and pure gluonic resonances, respec- 
tively. For the first case, as discussed above, we take 
T(R -> 77) = e 4 q (3a) 2 T tot (R)/2. Using the Particle 
Data Group (PDG) average values for the total width 
one gets r(/ (1500) -> 77) ~ 0.3 keV, r(/ (1710) -> 
77) ~ 0.4 keV and r(X(1835) -> 77) ~ 0.2 keV. The 
corresponding cross sections using these theoretical es- 
timates for the width are 3(158) /ib, 3.4(216) /ib and 
1.1(84) /ib at RHIC(LHC). If we are conservative, one 
can consider the experimental upper bounds for the two- 
photon widths of /o(1500) and /o(1710). This procedure 
gives an upper limit of the cross section for those reso- 
nances in peripheral collisions. The ALEPH experiment 



[3| studied the 77 production of those glueball candi- 
dates via their decay to ir + ir~ and the following limits 2 
were determined: r(/o(1500) — > 77) < 1.08 keV and 
r(/ (1710) -> 77) < 21.25 keV. Using those limits the 
corresponding cross sections are of order 0.95 mb (20 fib) 
for / (1500) and 11.5 mb (180 /ib) for / (1710) at LHC 
(RHIC) energies. We quote Ref. [1J] for a comparison 
of our results with a wide class of theoretical models and 
exotic QCD states in the meson production in photon- 
photon process. 

For a pure glueball resonance we follow [13}, adapted 
for the candidates considered here. Namely, assuming 
the /o resonances as states with J = L = then Eq. 
(54) of Ref. [l3j has been used, where we take the 
PDG values for the radiative J/ip decays in the fol- 
lowing channels: ip — > 7/o(1500) — > "/tttt and ip — > 
7/o(1710) -> -fKK. Assuming the X(1835) resonance 
to be a state with J = and L = S = 1 we relay on 
Eq. (35) of [l3j and use the PDG value for the de- 
cay channel ip — > 7^(1835) — > jir + ir~rj'. Putting all 
together, the estimates for the two-photon width for a 
pure glueball resonance are r(/o(1500) — > 77) ~ 0.77 eV, 
r(/ (1710) -^ 77) ~ 7.03 eV and r(X(1835) -> 77) ~ 
0.021 keV. Notice that for the X(1835) a larger width 
is predicted @, being of order 1.1 keV. The widths are 
about three orders of magnitude smaller that for pure qq 
states. Therefore, as the two-photon cross section scales 
as (2J+ l)r (R — > 77), Eq. ©, one can consider the ex- 
perimental feasibility of using peripheral heavy-ion colli- 
sions to determine the nature of the resonances discussed 
above. The values for the corresponding widths and cor- 
responding cross sections estimates are shown in Table 
I. 

Now, we address the Pomeron-Pomeron channel. In 
particular, we focus on the central diffraction (double 
Pomeron exchange, DPE) in nucleus-nucleus interac- 
tions. As a starting point we compute the DPE proton- 
proton cross section making use of the Bialas-Landshoff 
|15l fl6l | approach. We believe that this non-perturbative 
approach is a reasonable choice due to the light mass 
of glueballs candidates considered in present calcula- 
tion. For a perturbative QCD guided calculation we 
quote the recent work in Ref. [17j . where the exclu- 
sive scalar /o(1500) meson production is carefully inves- 
tigated. Here, we are particularly interested in the exclu- 
sive and central inclusive (central inelastic) DPE produc- 
tion of glueball states. In the exclusive DPE event the 
central object G is produced alone, separated from the 
outgoing hadrons by rapidity gaps, pp — > p + gap + G + 
gap+p. In the central inclusive DPE event an additional 
radiation accompanying the central object is allowed. In 
approach we are going to use, Pomeron exchange corre- 



1 For pure glueball resonance, G, the branching ratio is Br{G — >■ 
gg) ~ 1, whereas mixing states will give intermediate values of 
branching ratio. 



2 Here, we consider the ALEPH limits T(77 -> / (1500)) ■ 
Br (/o (1500) -» 7T+7I-) < 0.31 keV, 1X77 -> / (1710)) ■ 
Br (/o(1710) — > 7r+7r~) < 0.55 keV and taking the branching 
ratios 0.30 ± 0.07 and 0.026 ± 0.016 Q, respectively 



sponds to the exchange of a pair of non-perturbative glu- 
ons which takes place between a pair of colliding quarks. 
For DPE central inclusive G production we can neglect 
the additional gap spoiling effect, so-called Sudakov ef- 
fect. The scattering matrix is given by, 



\ a(t 2 )-l / \ a(ti)— 1 

M = M,,[j-\ (-J F(h)F(t 2 ) 



x exp {fi (fa + fa)) 5 gap (y/s) 



(4) 



Here Mq is the amplitude in the forward scattering limit 
(fa = fa =0). The standard Pomeron Regge trajectory 
is given by a (t) = 1 + e + a't with e « 0.08, a' = 0.25 
GcV~ 2 . The momenta of incoming (outgoing) protons 
are labeled by p\ and pi (fci and fe), whereas the glue- 
ball momentum is denoted by P. Thus, we can define the 
following quantities appearing in Eq. Q: s — (pi + P2) 2 , 

81 - (fcl + P) 2 , S2 = (fe + P) 2 , fa - (pi - fcl) 2 , 



fa 



(P2 



A:, 



The nucleon form-factor is given by 



F p (t) = exp(6i) with & = 2 GeV -2 . The phenomenolog- 
ical factor exp (/? (fa + fa)) with /3 = 1 GeV" 2 takes into 
account the effect of the momentum transfer dependence 
of the non-perturbative gluon propagator. The factor 
Sgap takes the gap survival effect into account i.e. the 
probability (S^ ) of the gaps not to be populated by 
secondaries produced in the soft rescattering. For our 
purpose here, we will consider Sg = 0.032 at y/s = 5.5 

TeV in nucleon- nucleon collisions 3 and S?. n = 0.15 at 

g a p 

y/s = 200 GeV (RHIC). In particular, for RHIC we have 
used an estimation using a simple one-channel eikonal 
model for the survival probability [lq . whereas for the 
LHC energy we follows Ref. [19| that considers a two- 
channel eikonal model that embodies pion-loop insertions 
in the pomeron trajectory, diffractive dissociation and 
rescattering effects. We quote Ref. [20| for a detailed 
comparison between the two approaches and further dis- 
cussions on model dependence of inputs and considera- 
tion of multi-channel calculations. 

Following the calculation presented in Ref. [16] we find 
.Mo for colliding hadrons, 



M = 32a 2 D 3 J d 2 Rp\vtp v 2 exp(-3K 2 /r 2 ), 



(5) 



where k is the transverse momentum carried by each of 
the three gluons. V^ v is the gg — > G J vertex depending 
on the polarization J of the G J glueball meson state. 
For the cases considered here, J — 0, one obtains the 
following result [lj| [2l[ : 



tiy&ti = 



2M 2 



A. 



(6) 



3 It is obtained using a parametric interpolation formula for the 
KMR survival probability factor |19(I in the form S'i = a/[b + 

ln(y/s/s )] with a = 0.126, b = -4.988 and s = 1 GeV 2 . This 
formula interpolates between CD survival probabilities of 4.5 % 
at Tevatron and 2.6 % at the LHC. 



where A is expressed by the mass Mq and the width 
T(gg — > G) of the glueball meson through the relation: 



A 2 =87rM G r(«7 ff ^G). 



(7) 



For obtaining the two-gluon decays widths the following 
relation is used, r (G -> gg) = Br (G -> gg)Y tot {G). 
At this point, some discussion is in order. The two-gluon 
width depends on the branching fraction of the resonance 
R to gluons, Br (G — > gg) and its knowledge would give 
quantitative information on the glueball content of a par- 
ticular resonance. As discussed before, it is a theoreti- 
cal expectation [12[ that Br (R(qq) —> gg) = 0(a 2 s ) ~ 
0.1 - 0.2 whereas Br (R{G) ->■ gg) ~ 1. Here, we will 
be conservative and assume the resonance to be a pure 
glueball. This fact translates into an upper bound for 
the exclusive DPE production as the cross section scales 
with r (R — >• gg). Following Ref. [22|], the two-gluon 
width can be computed from the resonance branching 
fraction in J/ip radiative decay, Br (ip — > 7 G) . For the 
candidates of interest here one obtains: 



Br (G(0++) -)■ gg) = 
Br(G(0- + )^gg) 



8tt(tt 2 - 9) Br[il> -)■ 7 G(0++)] M 2 



c R x\Hj(x)\ 2 T t 



^(tt 2 - 9)Br{iP ->• 7G(0"+)] M^ 



M 2 



c fl a;|Pj(a;)| 2 r tot 



M G ' 



where the function Hj((x) is determined in the non- 
relativistic quark model (NRQM) (see appendix of Ref. 
[22]) and cr is a numerical constant (Cr = 1, 2/3, 5/2 
for J PC = '", ++ , 2 ++ , respectively). The masses 
of J/ifj and of resonance are M^ and M G , respectively, 
and x = 1 — {Mq/MVj. Based on equations above, 
in Ref. [22| the following values for the branching 
fractions for scalar glueballs candidates are obtained: 
Br [/o(1500)] = 0.64 ± 0.11, Br [/o(1710)] = 0.52 ± 0.07. 
For the pseudoscalar X the situation is less clear due to 
small information on its decaying channels in radiative 
J/tp decays. The authors in [22[ have a prediction for 
77 resonance which gives Br [77(1410)] = 0.9 ± 0.2. As 
the branching fraction scales as 1/Mg in this theoreti- 
cal model, an educated guess for the X branching frac- 
tion would be Br [A(1835)] = (M^/Mx) ■ Br [r?(1410)] = 
0.69 ± 0.15. In the numerical calculations we set the 
limit case Br [X(1835)] = 1 and notice that the branch- 
ing would be about 30 % smaller. The values for T gg used 
in our calculations are summarized in Table II. A conse- 
quence on the small deviation for the branching fraction 
in pure qq and glueball resonance is the difficulty in test- 
ing their nature using the exclusive diffractive data. An 
option would be to obtain for instance the differential 
cross section on angular distributions and then compare 
the predictions for each composition (pure qq, mixing 
state and pure glueball). 

In addition, we use the parameters r = 1 GeV and 
D G 2 t = 30 GeV" 1 [3 where G is the scale of the 
process independent non-perturbative quark gluon cou- 
pling. An indirect determination of unknown parameter 



Glueball Candidate T 77 [eV] RHIC [nb] LHC [fib] 



h (1500) 


0.77 


14-9.3 


0.7-1.3 


/o (1710) 


7.03 


60-43 


3.8-8.6 


X(1835) 


0.021 


0.11-0.09 


0.01-0.02 



TABLE I: Cross sections for pure glueball candidates pro- 
duction through photon-photon fusion in electromagnetic 
nucleus-nucleus collisions at RHIC and LHC energies. 



a = G 2 /4:TT has been found in Ref. [23f using experi- 
mental data for central inclusive dijet production cross 
section at Tevatron. Namely, it has been found the con- 
straint S 2 ap (y/s = 2TeV)/a 2 l = 0.6, where 5g ap is the 
gap survival probability factor (absorption factor). Con- 



sidering the KMR [19j value 5| ap = 0.045 for CD pro- 
cesses at Tevatron energy, one obtains cxq — 0.274. 

The calculation presented above concerns to central 
inclusive process, where the QCD radiation accompany- 
ing the produced object is allowed. Therefore, in order 
to describe the exclusive processes where the central ob- 
ject is produced alone we include the Sudakov survival 
factor T(n,fi) [24] inside the loop integral over k. The 
Sudakov factor T(«, fi) is the survival probability that a 
gluon with transverse momentum k remains untouched 
in the evolution up to the hard scale fi = MqJ2. The 
function T(k,/j.) is given by 24]: 



T(k, fi) = exp 



M" a, | k 2 



dk 2 



1-8 



iP 2tt k 2 



zP 99 (z) + J2 P M 



dz , (8) 



where 8 = 



/(P 



P gg (z) and P qg {z) are the 



DGLAP spitting functions. In next section we will 
discuss the effect of introducing the Sudakov factor in 
the estimation of exclusive production in the Pomeron- 
Pomeron channel. 

In order to calculate the AA cross section the pro- 
cedure presented in Ref. [25] is considered, where the 
central diffraction and single diffraction cross sections 
in nucleus-nucleus collisions are computed using the so- 
called criterion C (we quote Ref. [251 ] for further de- 
tails). Using the profile function for two colliding nuclei, 
T AB = / d 2 bT A (b) T B (b - b), the final expression for CD 
cross section in AA collisions is given by |25| : 



°aa = A 2 J d 2 bT AA (b) exp [-A 2 a™ T AA {b)\ a 



CD 
PP " 



(9) 



where at™ and a 



CD _ 



gap 



- n ]: r Y{y/s) are the inelastic 



_ CD 

' pp "^^ ^ pp ^gap ' N PP 

and CD cross sections in proton-proton case, respectively. 
Using Woods-Saxon nuclear densities and considering 
the inelastic cross section a. 



(RHIC) energy, 



73 (49) mb for LHC 
5.5(0.2) TeV, we compute the 



/'/- 



Glueball r 99 [MeV] RHIC [mb] LHC [mb] 

/o(1500) 69.8 0.63 ±0.21 (inc.) 0.77 ±0.51 (inc.) 

0.40 ±0.14 (exc.) 0.50 ± 0.32 (exc.) 

/o(1710) 70.2 0.68 ±0.26 (inc.) 0.80 ±0.52 (inc.) 

0.41 ±0.16 (exc.) 0.49 ±0.31 (exc.) 

X(1835) 70.27 0.64 ± 0.24 (inc.) 0.77 ± 0.50 (inc.) 

0.38 ±0.14 (exc.) 0.45 ± 0.29 (exc.) 

TABLE II: Cross sections for inclusive (inc.) and exclusive 
(exc.) glueball production in double Pomeron exchange pro- 
cess for RHIC and LHC energies. 



CD cross section for nuclear collisions. The values for the 
inelastic cross section are obtained from DPMJET [26J J. 
where the scattering amplitude is parameterized using 
a tot , P and elastic slope (these parameters are taken as fit- 
ted by the PHOJET model [22]). We notice that for LHC 
energy the effective atomic number dependence is propor- 
tional to A 1 / 3 , which means that the nuclear CD cross 
section is only one order of magnitude larger than the 
nucleon-nucleon cross section. For completeness, we give 
the values of the DPE cross sections for the proton-proton 
case used in Eq. ©: a^(RHlC) = 0.170,0.180,0.168 
mb and o£ p D (LHC) = 0.80,0.85,0.83 mb for / (1500), 
/o(1710) and X(1835), respectively. 

In next section we compare the two production chan- 
nels and investigate the main theoretical uncertainties. 
We provide estimates of cross sections and event rates 
for both processes for RHIC and LHC energies at the 
heavy ion mode. 



III. RESULTS AND DISCUSSIONS 

In what follows the numerical results for the two- 
photon and Pomeron-Pomcron processes are presented 
and discussed. In Table I the cross sections for glueball 
production in photon-photon fusion at RHIC and LHC 
energies are shown. For RHIC we have considered the 
nominal center of mass energy of 200 GeV for gold-gold 
collisions and for LHC we take the planned nominal en- 
ergy of 5500 GeV in lead-lead collisions. The first value 
corresponds to the cross section obtained using a non- 
factorizable photon flux (Cahn- Jackson) (28[ and the sec- 
ond one refers to the factorizable flux as shown in Eqs. 
(PEJ). The deviation is sizable for RHIC and LHC. The 
cross sections are sufficiently large for experimental mea- 
surement. The event rates can be obtained using the 
beam luminosity [9j: for LHC one has £pbPb = 5 • 10 26 
cm s -1 , which produces the following number of events. 
One has 3.6 • 10 2 , 2 • 10 3 and 4 for / (1500), /o(lT10) and 
AT(1835), respectively, in the nominal LHC running time 
with ions of 10 6 s (one month). The event rates can 
be enhanced in a pPb mode, where the nominal beam 
luminosity is increased three order of magnitude com- 



pared to the PbPb mode. The present calculation can 
be compared to previous studies on glueball production 
in heavy ion collisions [29|, l3t| ■ In general, the numerical 
results are similar to those computations and the main 
deviation comes from the distinct estimates for the two- 
photon decays widths. A direct comparison can be done 
for the /o(1710) case, where in Ref. [30| one gets 48 
nb for RHIC and 2.3 /ib for LHC (using cut on impact 
parameter b > 2Ra and using T 77 ~ 4 eV [30J). 

For the convenience of phenomenologists we provide 
here a parameterization of the ultraperipheral AA cross 
section as a function of the resonance mass at the LHC 
energy. This makes simple the computation of event rates 
provided the specific meson state and its two-photon de- 
cay width. Using the Cahn- Jackson photon flux, we 
obtain in the interval 400 < M R < 4000 MeV the 
parametrization: 



ct M£ 



g upc (AA ->■ Rj + AA) ^ 

(2 J + 1) T{Rj -> 77) 1 + (M fl /4) ' 



(10) 



where ct = 4.9147 mb/GeV and /3 = -3.45335; T 77 and 
Mr are the decay width and the resonance mass in units 
of GeV, respectively. Several authors have argued for 
a low lying scalar glueball, with mass between 500 and 
1200 MeV [l|, [2|, depending on the proponents. The pa- 
rameterization above allows to obtain estimates starting 
from a modeling for the two-photon width. 

In Table II the results for Pomeron-Pomeron produc- 
tion of glueball is presented. The estimates are shown for 
the inclusive (inc.) and exclusive (exc.) double Pomeron 
exchange as discussed in previous section. Namely, for 
the inclusive production the Sudakov survival factor is 
not included (glueball is produced in association with 
Pomeron remnants) whereas for the exclusive case it is 
taken into account. In order to estimate the model de- 
pendence in the CD cross section, we have changed the 
soft Pomeron parameters in order to be consistent with 
the semi-hard Pomeron values considered in the DESY- 
HERA fits to diffractivc deep inelastic scattering (DDIS). 
For instance, taking FIT A of the HI Coll. [3l| parame- 
terization for the diffractive structure function F 2 one 
has e = 0.118, a' = 0.06 and b = 2.75 GeV~ 2 . Such a 
change enhances the cross section by a factor 3 for PbPb 
collisions at the LHC. In Table II, the cross sections are 
presented taking into account such a theoretical error 
band. Lower bound corresponds to soft Pomeron param- 
eters and upper bound stands for the semihard Pomeron 
ones. For RHIC energy, the Pomeron-Pomeron contribu- 
tion seems to be bigger than the photon-photon channel 
in a large extent. On the other hand, at the LHC they are 
competitive. However, the Pomeron contribution can be 
easily separated from photon channel by imposing a cut 
on the impact parameter of collision. After imposing this 
kinematic cut (b > 2Ra) the Pomeron contribution is re- 
duced as they are dominated by small impact parameter 
contributions. 

The present result is difficult to be compared directly 
to previous studies on Refs. [23, [3fJ. Those author did 



not include survival probability gap on their calculations 
and the theoretical approaches for Pomeron-Pomeron in- 
teraction are distinct. For instance, in Ref. [2{| the 
JPJP — >• G cross section is obtained using the Pomeron- 
quark coupling like a isoscalar photon, which allows to 
obtain the DPE cross section from the two-photon one. 
On the other hand, in Ref. [30j], only the inclusive double 
Pomeron production is considered. Following that study, 
we can perform a closer comparison. The cross sections 
are computed there with inelastic scattering effects us- 
ing the Glauber approximation (in Table 3 of Ref. [30] , 
see a aa elastic), which is similar to procedure presented 
here. After including gap survival probability factor one 
gets for the /o(1710) meson the values 1.23 (3.04) mb 
for RHIC (LHC), which is not so far from our results 
presented for inclusive production in Table II. 

Finally, it is important to discuss the uncertainties on 
the current calculations and the experimental feasibil- 
ity of detecting glueballs candidates. The main uncer- 
tainty here is the model dependence on obtaining the 
two-photon and the two-gluon widths for a pure glue- 
ball meson. For the two-photon width we considered 
a nonrelativistic gluon bound-state model of Ref. [L3| . 
which it could be a debatable issue and it is far from 
being optimal. There are more modern approaches as re- 
viewed in Ref. Q , but this is out of the scope of present 
work. For the two-gluon widths, we obtained them from 
the quarkonium width based on a non relativistic bound- 
state calculation [22[ . This type of matrix elements have 
been discussed in Refs. [32| giving rise to an effect of 
chiral suppression. We did not discuss the implication 
of those findings in present calculation. Concerning the 
experimental detection, the advantage of the exclusive 
processes discussed here is clear: glueballs are probably 
being produced with a high cross section in inelastic col- 
lisions (in pp or AA reactions) but when the multiplicity 
is high the combinatorial background is overwhelming. 
In exclusive production there is no combinatorial back- 
ground. In the ultraperipheral two-photon production of 
glueballs, the final state configuration is clear: nuclei re- 
main intact after collision and a double large rapidity gap 
between them is present (glueball is centrally produced 
with a low px transverse momenta spectrum). This type 
of measurement is already done at RHIC for photopro- 
duction of vector mesons and exclusive dilepton produc- 
tion with a signal identification well understood [33| . The 
situation for DPE glueball production is similar, with the 
Pt spectrum being broader than the processes initiated 
by two-photons. Thus, a transverse momentum cut (and 
also impact parameter of collision) could separate the two 
channels (for a review on these issues we quote Ref. [34]). 
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